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SUMMARY/OVERVIEW 


The  objective  of  the  present  program  was  to  study  the  structure  and  response  of  laminar 
premixed  and  nonpremixed  flames  with  emphases  on  effects  of  high  pressure,  flame/flow 
unsteadiness,  and  chemistry.  The  investigations  were  conducted  through  laser-based 
experimentation,  computational  simulation  with  detailed  chemistry  and  transport  descriptions, 
and  advanced  mathematical  analysis.  Phenomena  studied  include  the  development  of  intrinsic 
flameffont  pulsating  instability,  the  influence  of  aerodynamic  stretch  and  radiation  heat  loss  on  its 
development,  and  its  relation  to  extinction,  and  studies  related  to  the  development  of  detailed  and 
simplified  chemical  kinetic  mechanisms. 

Specifically,  the  results  show  that  in  the  absence  of  stretch  and  heat  loss,  the  onset  of 
pulsating  instability  is  described  well  by  the  Sivashinsky  criterion  of  Ze(Le-l)>l  1,  where  Ze  is 
the  Zeldovich  number  and  Le  the  Lewis  number.  The  presence  of  pulsating  instability  was  found 
to  narrow  the  flammability  limits.  The  onset  of  pulsating  instability  and  the  subsequent  flame 
extinction  occur  over  very  narrow  ranges  of  the  fuel  concentration  such  that,  while  its  effect  is 
substantial,  its  presence  cannot  be  readily  detected.  Furthermore,  it  is  shown  that  positive  stretch 
tends  to  promote  pulsating  instability  while  negative  stretch  retards  it,  which  is  complete  opposite 
to  the  trend  for  cellular  instability,  and  that  the  presence  of  radiation  heat  loss  further  narrows  the 
range  for  the  onset  of  instability  and  hence  the  flammability  limits. 

Studies  on  combustion  chemistry  have  been  proceeding  in  two  directions.  In  the  first, 
experimentation  was  performed  for  1,3 -butadiene  in  flow  reactors,  and  results  were  used  to 
compile  a  kinetic  mechanism  for  its  pyrolysis  and  oxidation.  Laminar  flame  speeds  also  were 
acquired  for  several  fuel  mixtures  at  atmospheric  pressure  and  methane/air  at  high  pressures.  A 
mixing  rule  was  developed,  emphasizing  the  importance  of  the  adiabatic  flame  temperature,  and 
the  adequacy  of  the  widely  used  methane  mechanism,  GRI-Mech,  was  scrutinized.  In  the  second 
effort,  mathematical  theory  and  computational  algorithms  were  developed  for  the  simplification 
of  detailed  mechanism  to  skeletal  mechanism  using  directed  relation  graph  tfieory  and  of  the 
skeletal  mechanism  to  reduced  mechanism  using  computational  singular  perturbation.  As  an 
illustration,  a  reduced  mechanism  was  derived  subsequently  for  ethylene,  demonstrating 
surprising  preservation  of  comprehensiveness  in  the  description  of  the  chemical  responses  under 
extensive  variations  of  the  thermodynamic  states  and  combustion  configurations  of  the  system. 

These  results  are  expected  to  be  useful  to  the  general  interests  of  AFOSR  in  the 
fundamental  emd  practical  issues  of  flame  dynamics  and  chemical  kinetics,  turbulent  combustion, 
soot  formation,  radiative  hear  transfer,  flame  extinction,  stabilization,  flammability,  and 
supersonic  combustion. 
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ACCOMPLISHMENTS 


The  research  program  has  two  major  focuses,  namely  the  unsteady  dynamics  of  flames  and  the 
chemical  kinetics  of  flames.  These  two  focuses  cover  a  substantial  arena  in  the  current  research 
activities  in  fundamental  combustion.  Good  progress  has  been  made  in  both  endeavors  during  the 
reporting  period.  Highlights  of  the  accomplishments  can  be  found  in  the  annual  reports  submitted 
to  the  program  director,  as  well  as  the  journal  papers  that  have  appeared  in  print  Thus  only  a 
very  brief  summary  of  these  works  are  mentioned  in  the  following. 

Unsteady  Flame  Dynamics 

Most  practical  combustion  situations  involve  turbulent  flames.  Under  most  situations  the 
chemical  reactions  take  place  in  narrow  regions,  constituting  the  so-called  laminar  flamelets. 
These  laminar  flamelets  are  subjected  to  the  fluctuating  shearing  action  of  the  turbulent  eddies 
and  under  severe  conditions  can  be  locally  extinguished.  Furthermore,  because  of  the  intrinsic 
disparity  in  the  difiusivities  of  the  heat  and  reactants,  the  flame  front  can  exhibit  intrinsic 
instability  in  the  forms  of  cellularity  and  pulsation  of  the  flame  surface.  The  presence  of  such 
instabilities  can  modify  the  burning  rate  and  possibly  even  affect  the  state  of  extinction.  The 
studies  on  effects  of  forced  oscillation  were  accomplished  primarily  during  the  previous  funding 
period;  consequently,  emphasis  was  placed  on  the  intrinsic  pulsating  instability  for  the  current 
reporting  period.  The  intention  is  eventually  to  couple  forced  and  intrinsic  instabilities  to  explore 
combustion  control  through  resonance.  The  following  is  a  synopsis  of  the  results  and  the 
understanding  gained  therein. 

Pulsatine  Instabilities  in  Premixed  Flames  [1.21 

It  has  been  well  established  that  premixed  flames  would  exhibit  pulsating  instability  for  Le>l 
mixtures  when  the  criterion  Ze(Le-l)>ll  is  satisfied,  where  Ze  is  the  Zeldovich  number  and  Le 
the  Lewis  number.  Since  Le  for  gas  mixtures  is  typically  close  to  1,  while  Ze  is  less  than  10,  it 
has  been  the  prevalent  belief  that  pulsating  instability  does  not  occur  for  gaseous  flames.  Rather, 
extensive  interest  has  been  placed  on  solid  flames  though  the  Self-Propagating  High-Temperature 
Synthesis  processes  because  Le»l  for  solids.  We  have,  however,  come  to  realize  that  because  of 
the  progressive  dominance  of  chain-terminating  reactions  relative  to  the  chain-branching 
reactions  near  the  states  of  extinction  and  flammability,  the  effective  activation  energy,  hence 
Zeldovich  number,  can  assume  substantially  larger  values,  leading  to  the  satisfaction  of  the  above 
criterion  with  the  potential  occurrence  of  pulsating  instability. 

In  this  study  the  onset  and  characteristics  of  pulsating  instability  in  unstretched,  planar, 
lean  heptane/air  [1]  and  rich  hydrogen/air  [2]  flames,  for  which  the  Le  is  greater  than  unity,  were 
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investigated  computationally.  The  occurrence  of  pulsating  instability  was  identified,  with  the 
pulsation  being  monochromatic  and  then  bi-chromatic  as  the  burning  becomes  weaker  with 
changing  equivalence  ratio.  In  the  presence  of  radiative  heat  loss  extinction  sets  in  shortly.  For  the 
heptane/air  flames,  the  transition  from  steady  to  pulsating  propagation  occurs  over  an  extremely 
narrow  range  in  the  fiiel  concentration  such  that  the  pulsating  extinction  mode  can  be  easily 
missed,  leading  to  inappropriate  interpretation  of  the  limits  of  extinction  based  on  steady-state 
extinction.  For  hydrogen/air  flames,  the  sharpness  of  the  chemical  kinetics  renders  the  transition 
between  steady  and  pulsating  propagation  quite  non-monotonic.  For  example,  at  elevated 
pressures  such  as  20  atmospheres,  increasing  the  fuel  concentration  can  lead  to  steady 
propagation,  pulsating  propagation,  and  steady  propagation  again  as  the  controlling  chemistry 
changes  from  that  of  H-O2  branching  to  HO2  branching. 

In  addition,  it  is  seen  that  the  pulsating  flames  are  quasi-steady  in  nature  in  that  the  period 
of  oscillation  is  larger  dian  the  characteristic  flame  time.  As  such,  the  unsteady  flame  cannot 
recover  once  the  instantaneous  flame  temperature  is  reduced  below  the  corresponding  steady- 
state  extinction  temperature. 

Pulsatine  Instabilities  in  Diffusion  Flames  [3,  41 

While  extensive  studies  have  been  performed  for  flamefront  instabilities  in  premixed  flames,  only 
recently  have  corresponding  studies  on  diffusion  flames  been  conducted.  The  reason  was  that 
diffusion  flame  sheets  are  inherently  stable.  However,  recent  analytical  studies  on  the  diffusion 
flame  structure,  allowing  for  finite-rate  kinetics,  have  identified  the  possible  onset  of  instability. 

Recognizing  the  potential  importance  of  such  phenomena,  we  have  eonducted  both 
computational  and  experimental  studies  with  real  fuel/oxidizer  systems,  aiming  to  scrutinize  the 
validity  of  the  analytical  findings.  The  studies  adopted  the  spherical  flame  as  the  configuration  for 
analysis  because  it  is  not  affected  by  aerodynamics.  It  was  determined  subsequently  that 
extinction  occurs  in  both  the  transport-induced  limit  of  low  mass  flow  rate  and  radiation-induced 
limit  of  high  mass  flow  rate.  Pulsating  instability  was  observed,  both  experimentally  and 
computationally,  near  both  limits.  The  pulsating  frequency  varies  from  a  few  Hz  for  methane 
flames  to  tens  of  Hz  for  hydrogen  flames.  Results  of  this  nature  are  expected  to  be  useful  in  the 
sensing  and  control  of  flames. 

Effects  of  Stretch  and  Heat  Loss  in  Pulsating  Instabilities  [5,  6,  7,  8,  91 

Effects  of  stretch  and  curvature  on  the  pulsating  instability  of  premixed  flames  have  been 
investigated  for  Le>l  mixtures  via  the  negatively  stretched  inwardly  propagating  spherical  flame 
and  the  positively  stretched  counterflow  flame.  Computational  and  asymptotic  analysis  results 
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show  that  pulsating  instability  is  suppressed  by  the  former  and  promoted  by  the  latter.  This 
behavior  is  completely  opposite  to  that  for  cellular  flames,  for  which  positive  stretch  tends  to 
suppress  instability,  while  negative  stretch  promotes  it.  In  particular,  for  a  given  rich  hydrogen/air 
mixture  whose  one-dimensional,  freely  propagating  flame  is  pulsatingly  unstable,  the  inwardly 
propagating  flame  initially  propagates  at  the  laminar  flame  speed  when  the  flame  radius  is  large. 
Oscillation  subsequently  develops  and  is  then  amplified,  damped,  and  eventually  suppressed  as 
the  flame  propagates  inward  and  the  magnitude  of  stretch  increases.  For  the  positively  stretched 
rich  hydrogen/air  counterflow  flame,  oscillation  is  initiated  at  an  equivalence  ratio  much  smaller 
than  the  one-dimensional  rich  threshold.  Furthermore,  the  critical  strain  rate  leading  to  pulsation 
is  smaller  than  the  corresponding  static  extinction  limit,  implying  that  the  flame  extinguishes  in 
the  pulsating  instead  of  the  steadily  propagating  mode  such  that  the  flammable  range  is  narrowed 
accordingly. 

The  computationally  determined  pulsation  limits  using  realistic  chemistry  and  transport  also 
were  found  to  agree  well  with  the  asymptotic  results  based  on  simplified  chemistry  and  transport, 
provided  that  the  global  activation  energies  and  Lewis  numbers  were  extracted  appropriately 
from  the  laminar  flame  responses. 

While  analytical  and  computational  results  show  that  pulsating  instability  occurs  for  Ze(Le- 
1)  =  11  »  1,  microgravity  experiments  conducted  elsewhere  showed  that  pulsating  instability 
also  occurred  for  sub-unity  Le  systems.  The  flammability  limits  also  were  narrowed.  A 
computational  study  allowing  for  radiation  heat  loss  then  demonstrated  that  this  loss  mechanism 
indeed  could  be  the  cause  for  the  narrowing  of  the  flammability  limits.  This  study  was  followed 
by  an  asymptotic  analysis  for  radiation-affected  stretched  flames.  The  experimental  and 
computational  results  were  captured. 

Flame  Propagation  and  Kinetics  [10, 11, 12] 

A  major  focus  of  this  program  is  to  acquire  experimental  data  on  flame  propagation  and  then  to 
use  these  data  to  validate  and  develop  kinetic  mechanisms.  We  have  been  engaged  in  studies 
involving  the  simpler  fuels,  including  hydrogen,  methane,  and  the  C2  hydrocarbons  ethane, 
ethylene,  and  acetylene.  Recently  we  also  have  performed  experiments  on  the  C3  hydrocarbons 
propane,  propene,  and  propyne.  For  the  present  period  we  have  moved  into  the  C4  group  of  fuels. 
The  first  one  we  attempted  was  1,3-butadiene  [10],  which  is  a  crucial  intermediate  in  the 
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oxidation  of  higher  hydrocarbons.  In  this  study  we  first  performed  flow  reactor  experiments  for 
its  pyrolysis  and  oxidation  over  the  temperature  range  of  1035-1185  K.  Using  these  data,  as  well 
as  those  from  the  literature,  a  kinetic  mechanism  for  1,3  butadiene  oxidation  was  developed. 

Considerable  effort  also  has  been  expanded  towards  acquiring  experimental  data  on  the 
laminar  flame  speed,  and  using  them  to  guide  the  development  of  kinetic  mechanisms.  Our  first 
effort  [11]  involved  the  determination  of  atmospheric  laminar  flame  speeds  of  mixtures  of  air 
with  ethylene,  «-butane,  toluene,  ethylene-n-butane,  ethylene-toluene,  and  n-butane-toluene. 
Binary  fuel  blends  with  1-to-l,  l-to-2,  and  2-to-l  molar  ratios  were  examined.  Experimentally, 
the  laminar  flame  speeds  were  determined  using  digital  particle  image  velocimetry  (DPIV). 
Since  the  use  of  DPIV  enables  the  mapping  of  the  two-dimensional  flow  field  ahead  of  the  flame, 
the  reference  speed  based  on  the  minimum  axial  velocity  point  as  well  as  the  imposed  strain  rate 
can  be  identified  simultaneously.  The  latter  now  can  be  determined  unambiguously  by  the  radial 
velocity  gradient  at  the  minimum  velocity  point.  By  systematically  varying  the  imposed  strain 
rate,  the  corresponding  laminar  flame  speed  was  obtained  through  nonlinear  extrapolation  to  zero 
strain  rate.  Computationally,  the  laminar  flame  ^eeds  were  simulated  for  all  single-component 
fuel/air  and  binary  fuel  blend/air  mixtures  with  a  detailed  kinetic  model.  Comparison  of 
experimental  and  computed  flame  speeds  shows  generally  good  agreement.  A  semi-empirical 
mixing  rule  was  developed.  The  mixing  rule,  which  requires  only  the  knowledge  of  the  flame 
speeds  and  flame  temperatures  of  the  individual  fuel  constituents,  is  shown  to  provide  accurate 
estimates  for  the  laminar  flame  speeds  of  binary  fuel  blends  under  the  conditions  tested. 

In  the  second  of  the  efforts  [12],  using  a  specially  designed  high-  and  constant-pressure 
combustion  chamber,  the  propagation  and  morphology  of  spark-ignited  expanding  spherical 
methane  flames  were  imaged  using  Schlieren  cinematography  and  a  high-speed  digital  camera. 
Stretch-ffee  laminar  flame  speeds  were  determined  subsequently  for  methane-air  flames  at 
pressures  up  to  20  atm.  and  methane-oxygen-helium  flames  at  pressures  up  to  60  atmospheres. 
Computational  simulation  using  GRI-Mech  3.0  showed  satisfactory  agreement  with  the 
experimental  data  up  to  20  atm.,  and  moderate  deviation  for  pressures  above  40  atm.  Markstein 
lengths,  global  activation  energies,  and  overall  reaction  orders  also  were  determined  as  functions 
of  pressure,  with  the  latter  two  parameters  exhibiting  non-monotonic  behavior  caused  by  the 
changeover  fi'om  H-O2  to  HO2  chemistry  similar  to  that  of  the  explosion  limits  of  homogeneous 
hydrogen-oxygen  mixtures. 
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Theories  and  Algorithms  for  Kinetic  Mechanism  Reduction  [13, 14, 15] 

When  the  present  research  program  was  initiated  a  few  years  ago,  the  stated  objective  was  the 
development  of  comprehensive  detailed  mechanisms  covering  hydrocarbons  all  the  way  from 
methane  to  the  C7-C8  range.  While  this  objective  is  still  very  much  the  main  focus,  an  additional 
goal  has  been  identified,  namely  the  development  of  mathematical  theories  and  computational 
algorithms  for  the  automatic  reduction  of  detailed  mechanisms  to  simpler  ones  without  loss  of 
comprehensiveness.  This  undertaking  was  motivated  by  two  considerations.  First,  detailed 
mechanisms  are  simply  too  large  for  integration  into  the  computation  of  combustion  flow  fields. 
For  example,  the  class  of  mechanisms  developed  by  Westbrook  and  co-workers  at  the  Lawrence 
Livermore  laboratories,  consisting  of  hundreds  of  species  and  thousands  of  reactions,  cannot  be 
used  even  for  the  calculation  of  the  one-dimensional  laminar  flame  structure  and  the  associated 
flame  speed.  Other  smaller  mechanisms,  such  as  GRI-Mech  for  methane  oxidation,  cannot  be 
used  for  turbulent  flame  calculations.  The  second  motivation  was  that  the  continuing  updating  of 
detailed  mechanisms  as  new  experimental  and  ah  initio  computational  data  become  available. 
Obviously  it  is  not  reasonable  to  expect  that  a  manual  reduction  needs  to  be  performed  whenever 
a  new  mechanism  appears. 

Our  foray  into  mechanism  reduction  was  initiated  a  few  years  ago  by  proposing  an 
augmented  reduced  mechanism  (ARM)  for  methane  oxidation.  This  mechanism  has  since  been 
used  extensively  by  many  turbulent  flame  and  engine  modelers,  with  satisfactory  results.  We  now 
have  extended  this  mechanism  to  include  NOx  formation  [13]. 

We  have  made  two  major  contributions  in  terms  of  theory  and  algorithm  development.  To 
understand  the  significance  of  these  contributions,  let  us  consider  the  general  procedure  adopted 
in  mechanism  reduction.  We  start  with  a  detailed  mechanism,  which  consists  of  many  species 
and  elementary  reactions.  The  calculation  time  can  be  estimated  roughly  as  being  proportional  to 
the  square  of  the  number  of  species.  This  detailed  mechanism  then  is  simplified  to  a  so-called 
skeletal  mechanism  by  eliminating  unimportant  species  and  reactions.  This  skeletal  mechanism 
then  is  simplified  further  to  a  reduced  mechanism  by  assuming  quasi-steadiness  for  some  species 
and  partial  equilibrium  for  some  reaction.  Thus,  while  the  skeletal  mechanism  still  consists  of 
only  elementary  reactions,  the  reactions  for  the  reduced  mechanism  are  lumped. 

The  crucial  criterion  in  mechanism  reduction  is  to  preserve  the  comprehensiveness  of  the 
chemical  phenomena  that  the  original  mechanism  is  supposed  to  describe.  The  reaction 
mechanisms  for  hydrocarbon  oxidation  are  intrinsically  very  complex,  being  highly  coupled  and 
nonlinear  in  terms  of  the  concentrations,  temperature,  and  pressure.  Another  well  known  example 
is  the  three  explosion  limits  of  hydrogen  oxidation.  For  a  given  temperature,  increasing  the 


8 


system  pressure  from  a  low  value  can  render  a  mixture  from  non-explosive,  to  explosive,  to  non¬ 
explosive  again,  and  to  explosive  again.  Furthermore,  the  mechanisms  are  expected  to  describe 
extensive  classes  of  combustion  phenomena:  from  atmospheric  combustion  in  furnace  and 
boilers  to  high-pressure  combustion  in  i.c.  engines;  from  low-temperature  ignition  to  high- 
temperature  burning;  from  lean  to  rich  reactions  as  the  mixture  equivalence  ratio  varies  due  to 
design  or  imperfect  mixing;  and  from  the  slow  smoldering  combustion  to  rapid  explosions. 

In  our  approach  we  first  generated  skeletal  mechanisms  from  a  detailed  mechanism  using 
directed  relation  graph  with  specified  accuracy  requirement  [15]  and  subsequently  generated 
reduced  mechanisms  from  the  skeletal  mechanisms  using  computational  singular  perturbation 
(CSP)  based  on  the  assumption  of  quasi-steady-state  species  [14].  Both  stajges  of  generation  are 
guided  by  the  performance  of  a  perfectly  stirred  reactor  for  high-temperature  chemistry  and  auto¬ 
ignition  delay  for  low-  to  moderately-high-temperature  chemistry.  The  demonstration  was 
performed  for  a  detailed  ethylene  oxidation  mechanism  consisting  of  70  species  and  463 
elementary  reactions,  resulting  in  a  specific  skeletal  mechanism  consisting  of  33  species  and  205 
elementary  reactions,  and  a  specific  reduced  mechanism  consisting  of  20  species  and  16  global 
reactions.  Calculations  for  laminar  flame  speeds  and  nonpremixed  counterflow  ignition  using 
either  the  skeletal  mechanism  or  the  reduced  mechanism  show  very  close  agreement  with  those 
obtained  by  using  the  detailed  mechanism  over  wide  parametric  ranges  of  pressure,  temperature 
and  equivalence  ratio. 

It  is  emphasized  that  Refs.  14  and  15  are  significant  contributions  to  the  theories  of 
mechanism  reduction. 

Review  Activities  [16, 17] 

Two  review  articles  were  written  on  the  two  major  activities  of  the  present  program.  Reference 
16  is  on  the  dynamics  of  laminar  flames.  It  updates  a  previous  review  article  on  the  same  subject 
presented  as  an  invited  paper  at  the  22”^  International  Combustion  Symposium  in  1988.  That 
article  has  enjoyed  extensive  citation,  by  both  the  laminar  and  turbulent  flame  communities.  The 
present  update  already  has  accumulated  an  impressive  list  of  citations.  The  second  article,  on  the 
development  of  detailed  and  reduced  mechanisms,  is  of  an  advocacy  nature,  bringing  attention  to 
the  importance  of  accurate  and  comprehensive  reaction  mechanisms  in  simulating  combustion 
and  propulsion  phenomena  to  the  aerospace  community. 

Writing  review  articles  is  a  major,  though  essential,  undertaking.  The  P.I.  is  pleased  to 
have  contributed  to  the  literature  in  this  manner. 
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SIGNIFICANT  INTERACTIONS 


1.  Contact:  The  Dow  chemical  Company;  Albert  Harvey;  (225)  353-4035 

Results:  Use  of  comprehensively  validated  reduced  chemical  kinetic  mechanism  allows 

realistic  description  of  methane  oxidation  chemistry 

Applications:  In-house  code  development  and  use. 

2.  Contact:  United  Technology  Research  Corp.;  M.D.  Colket;  (860)  610-7481 

Results:  Use  of  comprehensively  validated  reduced  chemical  kinetic  mechanism 

allows  realistic  description  of  meth2Uie  oxidation  chemistry 

Applications:  In-house  code  development  and  use 

3.  Contact:  Fluent,  Inc.;  Graham  Goldin;  (800)  445-4454 

Results:  Use  of  comprehensively  validated  reduced  chemical  kinetic  mechanism 

allows  realistic  description  of  methane  oxidation  with  NOx  formation 

Applications:  Commercially  available  code  for  engine  and  burner  simulation 
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